Quasi-Local Energy Flux of Spacetime Perturbation 
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A general expression for quasi-local energy flux for spacetime perturbation is derived from covari- 
ant Hamiltonian formulation using functional differentiability and symplectic structure invariance, 
which is independent of the choice of the canonical variables and the possible boundary terms one 
initially puts into the Lagrangian in the diffeomorphism invariant theories. The energy flux expres- 
sion depends on a displacement vector field and the 2-surface under consideration. We apply and 
test the expression in Vaidya spacetime. At null infinity the expression leads to the Bondi type 
energy flux obtained by Lindquist, Schwartz and Misner. On dynamical horizons with a particular 
choice of the displacement vector, it gives the area balance law obtained by Ashtekar and Krishnan. 
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I. INTRODUCTION 

It is well known that in General Relativity, because of 
equivalence principle, locally we can not detect gravity 
and therefore renders any unambiguous notion of local 
density for conserved quantities impossible. On the other 
hand, the notion of total conserved quantities for spatial 
and null infinity are well understood. Therefore, there 
has been hope that one should be able to find an appro- 
priate notion of quasi-local conserved quantities for finite 
spacetime domains [1]. 

A systematic way to study conserved quantities is 
through Hamiltonian. In a series of papers @, d, 0, 
a covariant Hamiltonian formalism was developed to ob- 
tain quasi-local energy-momentum expressions and found 
that the Hamiltonian boundary term also determines the 
boundary conditions for stationary spacetime. 

In general, for a dynamical gravitating system, the 
gravitational energy-momentum is not conserved, but in 
such case we would expect a meaningful total or quasi- 
local energy flux expression. In Q, covariant Noether 
charge (Hamiltonian) formulation was first used to iden- 
tify the Bondi energy flux. For trapping and dynamical 
horizons 0] , an energy flux expression was also obtained 
by Ashtekar and Krishnan Q . The formalism were then 
applied to trace the development of trapping horizons in 
black hole formation and evaporation Q. For general 
spacetime regions, definitions of quasi-local energy flux 
was also developed using covariant Hamiltonian formal- 
ism 0. The Bondi energy flux was derived but the en- 
ergy flux at dynamical horizons was not addressed there. 
Much discussions using marginally trapped surfaces and 
dynamical horizons were given in terms of the Vaidya 
spacetimes [Tfl ITTl fr3[. 
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In this paper, we shall construct explicitly a physically 
well defined energy flux formula for general dynamical 
gravitating system using the functional differentiability 
of the Hamiltonian and symplectic structure invariance 
of the theory which is independent of the choice of canon- 
ical variables and the boundary terms that one initially 
puts into the Lagrangian. We use spherical symmetric 
Vaidya spacetime as an example to test our energy flux 
expression 13]. The Bondi type energy flux at the null 
infinity and the energy flux at the dynamical horizons are 
also derived respectively in Schwarzschild and Painleve- 
Gullstrand coordinates. 

In section II, we derive our energy flux expression 
within covariant Hamiltonian formulation using func- 
tional differentiability and symplectic structure invari- 
ance. The important features of this expression are the 
independence of the choice of the canonical variables and 
the possible boundary terms that one initially puts into 
the Lagrangian for diffeomorphism invariant theories. In 
section III, we apply our expression to General Relativ- 
ity. Examples for Vaidya spacetime to obtain the Bondi 
type energy flux at null infinity are then given in the sec- 
tion IV. In there, the area dependence and implications 
of the first law is also observed. In section V, we calcu- 
late the energy flux for dynamical horizons and compare 
with the results obtained by Ashtekar and Krishnan Q 
for a particular choice of the displacement vector. We 
discuss the results in section VI. 



II. NOETHER CONSERVED QUANTITIES AND 
ENERGY FLUX 



In the following, we shall apply the Noether charge 
approach [E [E Ej, EE EE EE EE] for dynamical space- 
time using self consistence of functionally differentiable 
Hamiltonian and invariance of the corresponding sym- 
plectic structure to define a physical energy flux. Starting 
with a diffeomorphism invariant first order Lagrangian 4- 
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form C = C((j),p) — dip hp — A(ip,p), where <j> denotes 
an arbitrary collection of dynamical fields, and p being 
the corresponding conjugate momentum [l8|. The field 
equations, |§ = 0, ^ = are obtained by computing 
the first variation of the Lagrangian, 



S£ S£ 

SC = d(5<f> A p) + S<p A — + — A Sp 
0(f) op 



(1) 



where d(j) A p is the symplectic potential 3-form. For any 
diffeomorphism generated by a smooth vector field £, we 
can replace the variational derivative 5 by the Lie deriva- 
tive £^ 

£ C C = d{£f:<j) hp) + £/:(/) A — + — A £(: P , (2) 

Using the identity £^ = i^d+di^ and replacing di^C with 
£^C, one can define a conserved Noether current 3-form 
(or Hamiltonian 3-form Q) H(£) by 



H{i) = £ i (j>hp-i i L{ct>), 



(3) 



the Noether current 7i(£) is 



closed (dH{0 = -£^h^-^h£ s p ~ 0) when the field 



such that by equation 

equations are satisfied. Locally there exist a 2-form Q(£) 
(called the Noether charge) such that H(£) = dQ(£)- 
"field equation terms" or in general, 



no ^en+dQio. 



(4) 



Note that here are constraints including matter fields 
contribution. When integrated on a 3-space X, it gives a 
Hamiltonian 



h(o = [ no = [ en + 1 q, 

JS JS JdT. 



(5) 



therefore, Q can also be interpreted as the boundary term 
B = <f gs Q which defines the value of the Hamiltonian. 
We want to stress here, although H is called Hamilto- 
nian in the literature, it may not be functionally dif- 
ferentiable to define conserved quantities along the dis- 
placement vector field £ which generate diffeomorphism 
invariant transformations. Following the work of Regge 
and Teitelboim (l9j . we justify the functional differentia- 
bility of H as the total Hamiltonian by further varying 
H, 

SH(0= f SH(0= I S(£^Ap)-kS£(4>). (6) 
is Js 

Using equation ([1]) and dropping the field equation terms, 
we arrive at 



5H(£) = / 8{£^Ap)- £ s (54>hp)+ d> i s {S(j)hp) 
Jt, J os 



(£^4> A Sp - 5<j> A £%p) + f i s (5<t>hp), 

S JdX 



and the symplectic 3-form £1 is defined by 



^(5i,S 2 )= / 5i4> A 5 2 p - 6 2 (t> A Sip. 



(8) 



There can be two possibilities. If j> as ^(<50 A p) = 
then H(£) is automatically functionally differentiable, 
conserved along the vector field £, i.e. £^H(^) = 0. 
When one can find boundary conditions (i.e. see ref- 
erence (H) to give § ds i^(S(j) A p) = 5B(£), to modify 
H to H = H(£) - B(£) such that H is still function- 
ally differentiable, conserved along the vector field f, i.e. 
£^(0 = 0. The function H{£) and H{£) are called 
the functionally differentiable Hamiltonian conjugate to 

In general dynamical gravitating systems, when space- 
time is non-stationary, there does not exist boundary 
conditions to achieve a functionally differentiable Hamil- 
tonian to define conserved quantities. However, the re- 
placement of 5 by £^ in equation{7]) will still lead to the 
following flux expression, 



£^H(0 = I 

JdT, 



hp). 



(9) 



But this flux expression has the ambiguity on the choice 
of the canonical variables between cj) and p. Following [f| , 
a prescription was developed by insisting on the self con- 
sistency of functionally differentiable Hamiltonian and 
invariance of symplectic structure (symmetric under —<fi 
and p interchanges) [2(| • 

To proceed, we perform a second variation A on £^H 
in equation ©, and identify the energy flux as in the 
fallowings . We have, 

A£ ( H(£)=(£ i^(£^A(f) hp + £{<f) A Ap), (10) 
J os 

where A£^ = £ ^A is assumed. We subtract a Lie deriva- 
tive term on both sides, 

A£ i H(0-£i f «c(A<Mp) = f it(£t</>AAp-A(f>A£tp), 

JO'S JdT, 

(11) 

and observe that the last term has a symplectic structure 
and therefore allows us to define a quantity E(£) by 



(7) 



E(0=AJT(0-* i^hp). (12) 

JdY. 

Note that the above equation (fTI]) bears the same form 
as equation iff]), we arrive at the first main result of this 
work, 



F(0 = j£ e E(0 = i> i i (£ i <phAp-A(l)h£ i p). (13) 

From the surface integral form of F(£), non-conserving 
nature and only being functional differentiable to gen- 
erate the correct dynamical evolution on the surface 
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S = <9£, we can therefore interpret it being the total 
energy flux across some dynamical surface boundary S. 
To further confirm this energy flux interpretation, we ob- 
serve another important feature of this flux expression, 
namely, the symplectic structure invariant character will 
automatically gives the correct boundary terms, bound- 
ary conditions and a conserved Hamiltonian as discussed 
in Q in the stationary cases. 

Note that our general expression F(£) is precisely the 
sum of the two flux expressions of [4], 



F(f) = F D irichlot(0 +F Ncumann (£) 
— ^dynamic (£) H - -^constraint (£) , 



(14) 



where F Dirichlot (£), F Neumann (£), F dynamic (£) and 

^ constraint 

(£) are the flux expressions with certain 
variables being fixed on the boundary (see Appendix). 
For a general dynamic spacetime, we want to allow our 
symplectic variables to be completely dynamical without 
any variable being fixed on the boundary. An example is 
dynamical black holes which we shall discuss in section 
V. 

We would also like to point out that it is the functional 
differentiability and the symplectic structure invariance 
which frees the energy flux expression from the ambigui- 
ties in determining the boundary terms one initially put 
into the Lagrangian and therefore allows us to define the 
energy flux expression uniquely. 



III. ENERGY FLUX FROM SPACETIME 
PERTURBATIONS 

We now apply the prescription described in the previ- 
ous section to General Relativity, 



s = [ c = -±- [ ir fc A(*(tf a Atf fc )) + /; r , 



-> (15) 



where R ab = duj ab + ui a c A uj cb is the curvature 2-form 
constructed from the connection 1-form ui ab , *($ a At9 b ) = 
ie afc cc i$ c A and g = rj ab $ a <g) s i) b is the metric, where 
ifob = diag(— 1, 1, 1, 1) and d a is the orthonormal frame 
1-form field. It is important to note that one can add 
a boundary term to the above action, however, adding 
such a boundary term will not change the result using 
our symplectic structure invariant prescription described 
in the previous section. 

Denote 7 ]ab = (*(0 O A 6 )) = ±e abcd i9 c A uj ab being 
the spin connection that solves the equation of motion. 
The expression of Noether current 3-form in equation((5|) 
is given by 

H(0 = [ H(0 = [ dQ(0 = -^<f i^ ab A Vab , 

(16) 

where Q(£) is the Noether charge 2-form appears as a 
total derivative therefore can also be interpreted as the 



boundary term, here, we assume the field equations are 
satisfied. For stationary spacetime but rather generic 
situation {§ dTi i^4> Ap^O), we require 



5H = 



5(H - B) 
1 



16tt 
1 

167r 



8 <t> i^uj ab A r] ab 
as 



i^5u; ab Ari ab ) 



^w ab A Sr] ab + Suj ab A ii:r] ab 



= 0. 



(17) 



This can be satisfied by the boundary conditions on a 
bifurcate Killing horizons for stationary black holes [3, 
GJ, El such that £ £ tf° = 0, £ i oj ab = 0. In this case the 
energy flux is zero, £^H = and from equation (|16|) we 
can obtain the Noether charge for stationary case. 

For dynamical cases, following equation (|L3[) we obtain 
the corresponding energy flux formula using the previous 
described invariant symplectic structure prescription, 

= ^ri H ( £ i" ab A Af lab - Au ab A £ i 7 lab ) . 

(18) 

For small perturbation away from stationary spacetime, 
we can define Auj = u dynamical - ^stationary, and Al] = 
Vdynamicai - ^stationary ■ This is a consistent definition, 
as when in the stationary limit, the expression will give 
zero flux, correct boundary conditions and correct value 
of the Noether charges for stationary black holes. 



IV. VAIDYA SPACETIME EXAMPLE, BONDI 
TYPE ENERGY FLUX AND FIRST LAW 

As an example for calculations of our energy flux ex- 
pression, we consider the Vaidya spacetime which de- 
scribes a spherically symmetric collapse of null dust (ra- 
diation). The metric is given by 

ds 2 = ~e 2 ^dt 2 + e-^dr 2 + r 2 {d9 2 + sin 2 9d<j) 2 ), (19) 

where ip = ip(t, r) and e 2 ^ = 1 — 2m(i, r)/r. 

In this coordinate, the marginally trapped surfaces are 
given by r = 2m(t,r). For constant m(t,r), this is just 
the standard Schwarzschild metric. Now consider a per- 
turbation Am(t, r) away from the stationary solution, 



m(t, r) = niQ + Am(t, r), 



(20) 



because mo is a constant, this implies to' = d r (Am) and 
to = dt(Am). In terms of the orthonormal frames, the 
natural choice is, 



i9° = e+db, tf 1 = e'^dr, 
tf 2 = rd8, tf 3 = rsin 

with corresponding basis vectors: 



e = e-+dt, 

e 2 = -de, 
r 



ei = e <9 r , 
1 

e.3 = — —. 
rsm( 



(21) 



(22) 
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For £ = c\dt + C2d r , the Lie derivatives of $ a are 

dt, 
dr. 





= — -r c\m + C2 1 m - 


£ £ <& 1 


1 


c\m + C2 ( m' - 




£^ 2 


= c 2 d6, 




£^ 3 


= C2 sin 


9dcj). 



r 
rn 



(23) 



The spin-connection ui ab has the following nonvanishing 
terms: 



01 _ 




- l -(m'- 
r \ 


r / 


12 = 


-e^d0 = 






13 = 


-e* sin (9c 






23 = 


— cos Odd) 


= -c 32 . 





(24) 



The corresponding Lie derivatives of ui ab have the follow 
nonvanishing terms, 



1 



Ci / 4m 
' m -' 



■re 



2e 



rfr 



1 

£i 
r 

C2 



£2 

r 



. , Tn 
m 

r 



Amm! m 4mm \ 

7* y2g2^ J 



dt 



dr 



2m' ^ 2m 



p , , 12 



£^ 23 



1 

re^ 

1 

re tp 

P , ,31 



C\m + c 2 



cim + C2 



(m' 



m 
r 



m 
r 



dt 



dO 



= -£^lo 3 



(25) 



The perturbation of the orthonormal tetrad the spin- 
connection have the following forms respectively, 

Atf° = Urn*, 

re* 



Az? 1 = 

Ad 2 = 
Ai? 3 = 



1 



o, 

0. 



Am dr. 



(26) 



Ac; u 



/ 4mAr 



1 



re 



4i/' 



Am c?r 



, Am . , 
Am' ) dt 

r \ r 



-Auj 



10 



Auj 12 = —AmdO 



Aw 13 



Alj 



2:; 



1 

re* 
1 

re*~ 
= 



Am sin 
-Aw 32 . 



Auj , 
= -Ac 31 , 



(27) 



Many of the terms vanish in the energy flux F(£) , the 
remaining nonvanishing terms that will contribute are, 



— * i c (£^ 12 A A m2 + £(:uj 13 A Aijia) 

JOS 



Am r 



8tt 



=2^ 



Cim + c 2 



("•'-") 



(28) 



i^Aco 12 A £ 5 ?7i2 + Alj 13 A £ € ryi 3 ) 





Am 


( i m Y 


Am' 


Cl 




c\m + c 2 yn j 


- c 2 






r 



(29) 



and 



-!- <f ^(A W 01 A% 01 ) 

87r JOE 



f'2 



/9S 

. Am\ / 4m Am Am 
Am ~ — ) Cl -(^- + ^' f ^ 



F(0 



Finally, the total energy flux is 

2 Am \ /4mA/;/ 



(>2 



Am' 



Cl 



Am 



(30) 

C2 ■ 

(31) 



Taking u = t — r = const and t, r — > oo to approach 
the null infinity and dropping the term contains mAm 
which is of higher order in A, we arrive at the Bondi type 
energy flux 



= c 2 (cim' - c 2 m) = —d u m(u), 



(32) 



where we have put c\ = 1 and C2 = 1 by requiring 
that £ and A defines the same direction of mass changes 
for the consistency of interchanging A and £f. The 
same energy flux result, equation (|32p . was also obtained 
long ago by Lindquist, Schwartz and Misner [L3| using 
Landau-Lifshitz stree-energy pseudotensor. Such energy 
flux —d u m(u) has the interpretation as the luminosity of 
the star as seen by an observer at null infinity. 

It is interesting to notice that the nonvanishing term 
come only from the following equations 



no 



16tt 



i i [Auj 01 }A£ i (d 2 Ad 3 ) 



^-l fetf 1 )^ 12 A£<d 3 - Auj 13 A£rd 2 } 



16tt 



;/e 



(33) 
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where d 2 A$ 3 is the area element. This indicates the first 
law for general spacetime regions. For stationary space- 
times where £^('& 2 A t? 3 ) = 0, the energy flux vanishes. 
The appearance of the first law provides a nontrivial con- 
sistent check of our energy flux expression. 



V. ENERGY FLUX FROM DYNAMICAL 
HORIZONS: PAINLEVE-GULLSTRAND 
COORDINATES 

On dynamical black-hole horizons, r — 2m, which im- 
plies m! = 1/2, Am' = 1/2, m/r = 1/2. The first term of 
the above total energy flux (f3Tj) becomes 



no 



ClC 2 



1 - 



4Ar 



(34) 



Note that unlike in the Bondi type energy flux case which 
is defined at the null infinity 2-spherc boundary, here for 
dynamical horizons, £ is bounded by two cross sections, 
<9£ = dTi\ + 9£ 2 , with radius changes from r% to r 2 dy- 
namically because of the outgoing energy flux. However, 
the second term in (|31|) is singular on dynamical horizons 
where = (although m = 0) because of the coordinate 
singularity on the horizons. 

In order to study the energy flux in dynamical hori- 
zons, we make a coordinate transformation to the 
Painleve-Gullstrand time coordinate T which is related 
to the Schwarzschild coordinate t by, 



T = t + 4m 




(35) 



In terms of the Painleve-Gullstrand coordinates, the met- 
ric can be written as 



ds 2 = -dT 2 + (dr 



2 _ rr-> , ,.. : . l 2l7 2lIlll dT f + r 2 {de 2 +s - m 2 ed(j) 2y 



with the choice of the orthonormal frames, 



(36) 



tf> = dT, ^= d r + J 2 -^±dT, 



i9 2 = rde, i9 3 = r sin 



(37) 



and the corresponding basis vectors, 



2m(T,r) a 
e = d T - \ o r , ei = d r , 



e 2 = -d e , 
r 



&3 



(38) 



In order to compare results for dynamical horizons 
obtained by Ashtekar and Krishnan [8j previously, we 
choose an r dependent c 2 in the displacement vector, 



£ = ci<9t + c 2 {r)d r instead of a constant c 2 . We can 
now calculate the Lie derivatives of t9 a and obtain 



£^° = 0, 



V / 2mr 
£^ 2 = c 2 {r)d9, 
= 5 2 (r)sin 



(771 \ 1 
mf J dT + c f 2 (r)dr, 



(39) 



The corresponding spin-connection oj ab will have the fol- 
lowing nonvanishing terms: 



UJ 



01 



m ( 1 

-(m ) —f= 

r \v2mr 



dr 



1 



-dT 



,10 



UJ = — 



UJ = — 



2m 
r 

2m 



d6 = -uj 2a , 



sin 



,30 



UJ 



12 



UJ 



13 



uJ 



23 



-de -- 

- sin 

- cos 



,21 



,31 



-U> 



32 



(40) 



and the Lie derivatives of oj ab have the follow nonvanish- 
ing terms: 



£, 



LU 01 = 



, mm m , , 

( " m+ ^T + 27 )dr 



y/2mr 

C2(r) , „ , (m 



(-m" + 



77? 
7" 



3m 
2^2 



V2mr v 2m 

— ==( m)dr + ci(— )dT 

V2mr r r z r 

, . , ml m m" . ,„ 
- C2 ( r) (2--2--— )dT 



7° , ,03 



£ f uj LZ 
£ f uj 13 
£ f uj 23 



1 



V2mr 
7 1 , ,20 

1 



V2T 



cim + c 2 



cim + c 2 



(r) m 

( \ ( * m 
(r)[m 



de 



= -A 
= o = 
= o = 

= = 



,30 



-£^ , 

-£^uj 32 . 



(41) 



As before, we express the perturbation of mass as, 

m(T,r) = m Q + Am(T,r), (42) 



A 



Am 



1 -Am 
•>/m 2m 1 /m' 



(43) 



(44) 
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and obtain the following perturbation of the orthonor- 
mal tetrad and spin-connection respectively, 



Ai9° = 0, 
Atf 1 



1 



V / 2mr 
Atf 2 = 0, 
Atf 3 = 0, 



(Am)dT, 



(45) 



Alo 01 = -Aw 10 

, . , Am, / 1 1 , 

= -(Am' ) - dr + -dT 

r \\j2mr r 

. . m. Am 

+(m' ) ; dr, 

r 2my/2mr 

Auj 02 = ^={Am)de = -Alo 20 , 

Aw 03 = -L=(Am)sin0e^ = -Aw 30 , 

\f2mr 

Au 12 = = -Aw 21 , 

Aw 13 = = -Au; 31 , 

Alo 23 = = -Alu 32 . 



(46) 



Put all the above results into the energy flux expression 
F(£), we obtain the following nonvanishing contributions, 
which includes: 

— / i 6 (£(:uj 02 A A?7o2 + £^ 03 A A7703) 

= -l Am 1 5l52(r)m ' _ £i£aW^| ; (47) 

2 V m m r / ' 



--!- / i 5 (Aw 01 A £^01) 
= c 2 (r) 



+c 2 (r) 



(-'-7) ( £2 «^\/^) 



(48) 



and 



c 2 (r)Am / Fhn _ , . \ 

+ ^2=r( v Cl Vv +c ^J 
+^(-'-7) (^/^) 



i A m5 2 (r)5 2 (r)^- 



(50) 



Note that an expression on the dynamical horizon was 
obtained by Ashtekar and Krishnan [8j previously, 



t-t) =L T - bT ° ( ^ v 

Af r {|,r| 2 + 2K| 2 } d'V. 



I AH 
1 



16tt 



AH 



(51) 



At the dynamical horizon, r = 2m, therefore m' = 
1/2, Am' = 1/2, m/r = 1/2, our flux expression (|50|) 
reduced to 



F(0 



, ^ , , , 4 Am, 
c 2 (r)(c 1+ c 2 (r))(l ) 



-- [Amc 2 (r)c 2 (r)]^ . 



(52) 



Therefore, on the dynamical horizon when we choose 
ci = and c 2 (r) = y^r for the displacement vector £, 
contribution from the second term vanishes and our re- 
sult reduced to an "area balance law" , 



(53) 



which agrees with the Ashtekar-Krishnan energy flux for- 
mula (|5"Tj) . with the shear \a\ and the twist |£| becomes 
zero in Vaidya spacetime. 

Note that similar to the previous section, in this ex- 
ample the nonvanishing term come from 

F(0 



16tt 



z c [Aw ul ] A£ c (z9 2 Atf 3 ) 



~ <f i e (Acj 02 A £^ Q2 ) -^-i i^A^ 3 A % /03 ) 
c 2 (r)Am 



V2~r 



5l V ~r" + 



— I {^^[Alo 02 A£ s $ 3 - Auj 03 A£^ 2 } 

1 / . r * m Am 

— * z e [-Aw 01 + 

D7T J as ry/2mr 



1,1 ' Am ^]A£ $ (^A#), 



1 / c 2 m ci 62 (r)m' cic 2 (r) 
—Am 1 

2 \ m m r 

-iAm5 2 (r)5 2 (r)X^. 



Finally, the total energy flux becomes, 



(49) 



(54) 

where t9 2 A t? 3 is the area element and again hints at 
an area dependent first law. For stationary spacetimes 
where £^{'d 2 A i9 3 ) = 0, the energy flux vanishes. 



VI. DISCUSSIONS 

A general expression for quasi-local energy flux expres- 
sion is derived from covariant Hamiltonian formulation 
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using functionally differentiability and symplectic struc- 
ture invariance, which is coordinate independent. The 
energy flux expression is given by the boundary term. 
The benefits of using symplectic structure invariance is 
to avoid the ambiguity in choosing the correct bound- 
ary terms in the Lagrangian that one begins with. This 
was a core problem for many other formalisms for diffeo- 
morphism invariant theories in the literatures. Another 
important features of this expression are the indepen- 
dence of the choice of the canonical variables. The ex- 
pression F(£) depends on the vector field £ and the choice 
of the boundary surfaces which depends on the type of 
the physics under investigation. 

For the boundary surface taking to be the null infinity, 
the expression leads to the Bondi type energy flux ob- 
tained by Lindquist, Schwartz and Misner [131 ] for Vaidya 
spacetime, where the energy flux has the interpretation as 
the luminosity of the star as seen by an observer at infin- 
ity. If the boundary surface is taken to be the dynamical 
horizons, the expression gives rise to the energy flux ob- 
tained by Ashtekar and Krishnan for Vaidya spacetime. 

Note that for Painleve-Gullstrand coordinates our ex- 
pression fails to give the correct Bondi type energy flux 
at null infinity because of the 1/y/r fall off of the metric. 
Also for Schwarzschild coordinates the expression fails on 
the interesting dynamical horizon. This only indicates a 
trivial fact that one requires to use at least two coordinate 
patches to cover the whole spacetime range of interest ex- 
cept for the uninteresting flat spacetime. However, this is 
a good news to our energy flux expression which derived 
from a coordinate independent covariant formalism and 
detects the limits of the applicability of the coordinate 
system employed in the calculations. 

Another interesting observation is that the expression 
gives an area dependence which hints at the first law for 
general Vaidya spacetime. 

Appendix 

In comparison to the flux expressions in [4], our flux 
expression (IT5]) is simply the sum of the Dirichlet and 
Neumann boundary flux expressions given in [4]. In the 



special case when <f> is a fixed variable, A<fi = 0, our flux 
expression (Ti~3"|) reduced to the "Dirichlet boundary flux 
expression" (equation (30) of [4]), 

F D irichlet(0 = f k(£^AAp). (55) 
JdY, 

If p is a fixed variable, Ap = 0, our flux expression 
(I13p reduced to the "Neumann boundary flux expression" 
(equation (31) of [4]), 

F Neumann (0 = f if(-A<M £ 6 p). (56) 

JdY, 

Similarly, when the spatial projections of the variables 4> 
and p are fixed, we obtain the "dynamic boundary flux 
expression" (equation (32) of [4]), 

Fdynamic^) = / (tftf A ^Ap - i f A^ A (57) 
JdT, 

where ? = (— l)' for <fi being an /-form field. If the time 
projections of the variables (f> and p are fixed, we obtain 
the "constraint boundary flux expression" (equation (33) 
of [4]), 

F C onstraint(0 = f (^^^^ A Ap - <rA ? A i ? £ ? p). (58) 

J as 

Our flux expression (|13p is the general symplectic invari- 
ant expression without any variable being fixed. 
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